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Crystal structures of CLi 4 compounds are explored through ab initio evolutionary method¬ 
ology. Phase transition from metal to semimetal and semiconductor, and eventually to in¬ 
sulator with increasing pressure are revealed under pressure. Pressure-induced evolution of 
anti-metallization has been described quantitatively by Fermi Surface Filling Parameter and 
electron energy band gap using ab initio GW calculations. Anti-metallization is attributed to 
the hybrid valence electrons and their repulsion by core electrons into the lattice interstices. 
Very weak electron-phonon coupling interactions are found in the metallic phases, resulting 
in very low superconducting temperature. 

Hydrogen is the most abundant element in the universe. Though, at the top of the alkali 
metal column in the periodic table, hydrogen is not an alkali metal and is an insulator in ordinary 
conditions. The pressure-induced transformation is expected to be accompanied by a decrease in 
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the band gap and eventual formation of metallic hydrogen \ attracting a lot of interest. Lithium 
has a single valence electron similar to hydrogen, but lithium is a good conductor of heat and metal 
under ambient conditions. Unlike the metallization of hydrogen under high pressure, lithium can 
transform from a nearly free-electron metallic solid to an insulating one, i.e. pressure-induced anti¬ 
metallization 2 . As a contrasting counterpart, lithium can provide meaningful referential aspect for 
understanding of hydrogen at pressures. 

Recently hydrogen-rich materials have been investigated extensively and profoundly. Alkali 
metals, e.g. lithium and potassium, together with hydrogen can be compressed into alkali high- 
hydride alloys®, and metallize at pressures much lower than that required for pure hydrogen due 
to the ’’chemical precompression”. A special mention has been made on the group-IV hydrides, 
calculations and experiments on compressed silane 5 7 , disilane®, germane^, etc. show the pos¬ 
sibility of metallization at moderate pressures. These studies are believed to be useful to realize 
metallic hydrogen, which has been described as ’’the holy grail of high-pressure physics” ®. Ac¬ 
cordingly, the studies of group-IV lithium compounds should benefit to the comprehension of the 
metallization of group-IV hydrides which are very important for understanding and researching 
metallic hydrogen. 

Lithium is well-known for its ability to form heterogeneous clusters with various elements. 
Since the experimental discovery and theoretical verification of ’’hyperlithiated” bonding, which 
involves formal violations of the octet rule in doped Li clusters there have been lots of the¬ 
oretical and experimental investigations with the purpose of insight into the specific structural 
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and electronic properties of these systems. Such systems include CLi„ EES etc. CLi 6 was first 
predicted 13 and later verified experimentally ^ to be stable with respect to dissociation into the 
stoichiometric CLi 4 and Li 2 molecules. To enhance the H 2 binding and attain high storage ca¬ 
pacity, some nanostructures have been functionalized with polylithiated molecules, e.g. CLi 4 for 
hydrogen storage. The stabilities of the designed functional materials for H 2 storage have been 
verified by means of molecular dynamics simulations^ 15 [ 

In view of interests mentioned above, we have explored the crystal structures of lithium- 
rich compound CLi 4 in a wide pressure range from ambient pressure to 300 GPa. Contrary to 
the group-IV hydrides, a reverse pressure-induced metallization is found in the CLi 4 system under 
high pressure. We have found that the phase transitions from metallic PI, P-1, P2 1; C2/m, to 
semimetal and semiconductor of R-3m, and eventually to insulated Fddd with increasing pressure 
quantitatively by Fermi Surface Filling Parameter and electron energy band gaps using ab initio 
GW calculations. Pressure-induced anti-metallization is found not only in the CLi 4 system, but 
also in the phase of R-3m. Unlike group-IV hydrides, the electron-phonon coupling interactions 
are fairly weak and the superconducting temperatures are very low under high pressure. Finally, 
we attribute anti-metallization to the hybrid valence electrons and their repulsion by core electrons 
into the lattice interstices. 

The enthalpy-pressure curves of CLi 4 show the possible decomposition to the elements and 
reported compounds, see Figure 1, and the other possible lithium-rich carbides CLi n (n=2~6) and 
their stability can be found in supplementary material (see Extended Data Figure 1). Obviously, 
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the compound of CLi 4 can be stable at least above 1 bar and below 300 GPa. We carry out un¬ 
constrained searches at fixed pressures of lbar, and 5, 10, 50, 100, 200, 300 GPa using one, two, 
three, four, five, and six CLi 4 formula units if.u.) per cell, and have found many different relaxed 
structures. A special one with twelve units, i.e. sixty atoms per cell, is also performed at the low 
pressure of lbar and 5 GPa, respectively. Figure 1 shows our calculated enthalpies of the candi¬ 
date structures including USPEX and XY 4 -type ones, together with the decomposition (C+4Li and 
(C 2 Li 2 )/2+3Li) enthalpies (black and red lines, respectively). Before about 3 GPa, CLi 4 is unstable 
and decomposes into a C 2 Li 2 and Li mixture. After about 3 GPa, CLi 4 compounds become sta¬ 
ble energetically over the decomposition enthalpies lines, as shown in Figure 1. Several different 
structures with space group of PI are found at about 5 GPa with very close enthalpies. So, CLi 4 
compounds can be considered as amorphous at this low pressure range. To simplify the structure 
of amorphous phase, a small cell with thirty atoms is chosen as an example for the PI structures 
shown in supplementary material (see Extended Data Figure 3. a), and will be discussed below. 
Details of two PI structures as shown in Extended Data Figure 3 are provided in the supplementary 
material (see Extended Data Table 3 and 4). The P-1 phase becomes stable above 7 GPa, which 
is obviously energetically favored over the others until the pressure close to 18 GPa. At 18 GPa, 
P-1 phase starts to transform into P2| CLi 4 . Then P2, phase has the obviously lowest enthalpy 
among the structures until 40 GPa. After 40 GPa, P2i phase transforms into C2/m phase. At 60 
GPa, a high symmetry phase of R-3m appears and is stable until 128 GPa. The magnified images of 
enthalpy curves in this pressure range can be found in Extended Data Figure 2. A pressure-induced 
symmetrization processes from low symmetry phases to the high symmetry one with increasing 
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pressure until 128 GPa is demonstrated. After 128 GPa, another high symmetry phase of Fddd 
appears with 32 symmetry operators, which is slightly less than the 36 operators of R-3m. Obvi¬ 
ously, the five phases, i.e., P-1, P2i, C2/m, R-3m and Fddd, are the most favored structures of CLi 4 
compounds under high pressures. 

The favored structures of CLi 4 obtained in this work are shown in Figure 2, and the corre¬ 
sponding parameters at respective pressures are listed in Extended Data Table 1. At low pressure 
range, the triclinic phase of P-1 consists of two f.u. of CLi 4 (i.e., C 2 Li 8 ) with ten unequivalent 
atoms in the primitive cell, as shown in Figure 2 a. Each unequivalent atom occupies the crystal¬ 
lographic 2 i position and the site symmetry is 1. The P2 4 phase consists of two f.u. of CLi 4 (i.e., 
C 2 Li 8 ) with ten unequivalent atoms in a monoclinic crystal lattice, as shown in Figure 2 b. The 
eight Li atoms and two C atoms occupy the crystallographic 2 a position with the site symmetry of 
1. From the Figure 2 c, we can see that C2/m CLi 4 is another monoclinic phase. Consisting of two 
f.u. of CLi 4 (i.e., C 2 Li 8 ) in the crystal lattice, there are three unequivalent atoms in the conventional 
cell. Two Li atoms occupy the crystallographic 4 i position with m site symmetry, and one C atom 
locates at 2c position with 2/m site symmetry. Lattice parameters of this monoclinic phase at 40 
GPa are shown in Extended Data Table 1. The R-3m phase contains thr eef.u. of CLi 4 (i.e., C 3 Lii 2 ) 
in a hexagonal crystal lattice, as shown in Figure 2 d. The corresponding lattice parameters of 
this hexagonal phase at 60 GPa are listed in Extended Data Table 1. There are three unequivalent 
atoms, including two Li atoms and one C atom occupying the crystallographic 6c and 3b positions 
with 3m and -3m site symmetry, respectively. The Fddd phase consists of eight f.u. of CLi 4 (i.e., 
C 8 Li 32 ) in an orthorhombic crystal lattice, as shown in Figure 2 e. The corresponding lattice pa- 
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rameters of this phase at 150 GPa are listed in Extended Data Table 1. There are two unequivalent 
atoms in the conventional cell. The Li atom occupies the crystallographic 32 h position with the 
site symmetry of 1, and C atom locates at the 8 a site with the site symmetry of 222, respectively. 

After the discussion on the viewpoint of thermodynamics above, the studies focus on me¬ 
chanics and the lattice dynamics to test the reasonableness of propositional structures. The me¬ 
chanical stability provides a useful understanding for the structure of crystals. The strain energy of 
a crystal must be positive against any homogeneous elastic deformations, i.e., the matrix of elastic 
constants C l3 must be positive definite 29 . The elastic constants have been calculated and listed in 
Extended Data Table 5. Obviously, the elastic constants of the structures satisfy the mechanical sta¬ 
bility criteria 29 3 ^, indicating that the five structures are mechanically stable. The charts of phonon 
band structure and the projected DOS of P-1, P2i, C2/m, R-3m, and Fddd phases calculated at 
selected pressures are shown in Extended Data Figure 4. The absence of imaginary frequency 
modes indicates that these structures are stable dynamically. A direct result from PVDOS of these 
structures shows that the lighter Li atoms contribute more in frequency vibrational modes than the 
heavier C atoms. 

It is known that ab initio DFT calculation usually gives incorrect estimate of fundamental 
energy band gap of materials Recently, ab initio GW calculations show good improvement in 
predicting energy band gap. The GW correction can yield a reliable band structure of candidates to 
ensure that they are indeed metallic under pressure. In view of this, it is worth to perform electronic 
band structure calculations at the level of ab initio GW method. From the insets e and f of 3DFS of 
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R-3m at 60 GPa in Figure 3, we can clearly find the difference between ab initio DFT calculation 
and ab initio GW method. 

To investigate the degree of metallization quantitatively in metal CLi 4 under pressure, Fermi 
Surface (FS) Filling Parameter ( is used, where C is described as ( = jj J2kn £ kn — £f) °c c!/^ ' 
Obviously, this parameter can be deduced by the nesting function £ (Q) 32 . Here A(0) is the FS 
density of states, and C^fe^, accounts for the intersection point between envelope and Fermi line in 
the chart of electron density of state of metal. This parameter can be used to describe free electrons 
ratio in metal. Obviously, £ is the maximum value of £ (() ') at Gamma point and means the degree 
of metallization, as shown in inset g of Figure 3. 

A remarkable property of CLi 4 is the existence of phase transitions from metal to insulator 
with increasing pressures. At the lower pressure range below 70 GP, PI, P-1, P2 4 , C2m, and R-3m 
phases are all metallic, and the degree of metallization declines gradually with increasing pressure 
except a slight pressure-induced metallization from phase P-1 to P2 4 , as described in Figure 3 
(left part). At about 70 GPa, the metallic character of R-3m CLi 4 disappears, and the non-metallic 
character of energy band gap begins to appear. Above 70 GPa, the energy band gap of R-3m phase 
grows bigger and bigger until reaching a value of 2.43 eV at 128 GPa with GW corrections, as 
shown in Figure 3 (right part). So, the R-3m CLi 4 is a transitional phase from metallic property to 
non-metallic one. After 128 GPa, insulated Fddd phase appears with a band gap value of 3.26 eV 
at 128 GPa. With further increase of pressure, the band gap has a slight elevation until 300 GPa. 
So, the decrease of metallization and increase of non-metallic character are revealed in compressed 
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lithium-rich compound CLi 4 system, as described quantitatively by the FS Filling Parameter ( and 
electron energy band gaps. 

Subsequently, calculations about electronic structures have been performed to reveal the 
emergence of anti-metallization in the dense lithium-rich compound CLi 4 . Hybridizations of va¬ 
lence electrons (Li-s, Li-p, and C-2p) have been discovered under pressure, see PDOS in Figure 
4 (left part). At lower pressure range, the interaction between valence electrons and ionic cores is 
weak, and the nonlocal electrons near Fermi energy account for the metallic property, as displayed 
in Figure 4 a (middle and right). From 3D-ELF and 2D-ELF in Figure 4 a, we can see that the 
conductivity comes from connected regions where the value of ELF is about 0.5. But under higher 
pressure, volume contraction in crystal lattice increasingly shortens interatomic distances, and hy¬ 
brid electrons are repulsed by the core electrons into the lattice interstices, see Figure 4 b (middle 
and right). From 3D-ELF and 2D-ELF in Figure 4 b, we can see that nearly free electrons with the 
ELF value of 0.5 become totally localized and are confined to disconnected areas in interstitial po¬ 
sitions of atoms. As a direct result of charge distributions under pressure, charge transfer has been 
observed, see Extended Data Table 6. With the increased overlap of the valence wavefunctions 
with pressure, the s-state electrons transfer to the p-state electrons in element Li by compression. 
Similar charge transfers from s-state to p-state have been reported in pure lithium under high pres¬ 
sure 33 . Meanwhile, charge transfers from C-2s to C-2p, and partially to Li-2p have been observed 
by the integration of PDOS. So, high pressure induces the decreased electronegativity of element 
C and increased the one of element Li. 
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Finally, superconductivities are investigated under the framework of BCS theory, as show 


in Figure 3 (left part). The electron-phonon coupling parameter (EPC A), the logarithmic average 
phonon frequency (cei og ), and the Eliashberg phonon spectral function (a 2 F(u)) are calculated at 
high pressures. The resulting EPC A for PI, P-1, P2 1? C2/m, and R-3m phases at 5, 15, 30, 40, 
50, 60 GPa are 0.533, 0.398, 0.425, 0.364, 0.276, and 0.212, respectively, indicating that the EPC 
is very weak. The superconducting critical temperature can be estimated from the Alien-Dynes 


modified McMillan equation T c = ex P 


1.04(1+A) 


which has been found to be highly 


A—0.62A) 

accurate for many materials with A < 1.5. The cj] og can be calculated directly from the phonon 
spectrum, and the values for the metallic phases mentioned above, are 305.75, 416.16, 400.64, 
405.86, 484.95, and 560.85 K at selected pressures, respectively. The Coulomb pseudopotential 
/i* is often taken as 0.1 for most metals. Using n* of 0.1, the estimated T c are 4.725, 1.691, 2.344, 
0.949, 0.094, and 0.001 K, respectively, revealing a similar varied tendency with anti-metallization 
under high pressure. 


The dens lithium-rich compound CLi 4 demonstrates an unexpected process of anti-metallization, 
which is completely contrary to their counterpart of group-IV hydrides. Pressure-induced anti¬ 
metallization such as reported here also is expected in pure elements and compounds, such as al¬ 
kali metal 2 — and some alkali metal compounds, where the interaction between valence electrons 
and ionic cores is very strong. During metallizing or anti-metallizing in metallic states, Fermi 
Surface Filling Parameter (FSFP) will be a valuable parameter to quantify the evolution of the free 
electrons, as described here. 
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METHODS SUMMARY 


Ab initio evolutionary algorithm designed to search for the structure with the lowest energy at a 
given pressure has been implemented in USPEX code 16 . The underlying ab initio structure relax¬ 
ations have been performed using density functional theory within the Perdew-Burke-Emzerhof 
(PBE) parameterization of the generalized gradient approximation (GGA) as implemented in 
the Vienna ab initio simulation package VASP code ' 18 . The all-electron projector-augmented wave 
(PAW) method *^is adopted with the PAW potentials taken from the VASP library, where Is2s2p 
and 2s2p are treated as valence electron configurations for Li and C atoms, respectively. A plane- 
wave basis set with an energy cutoff of 1000 eV is used and gives well converged total energies. 

The plane-wave pseudopotential method within the PBE-GGA, through the Quantum-ESPRESSO 
(QE) package 2() ' is employed to study the electronic properties, lattice dynamics, and electron- 
phonon coupling (EPC). In the QE code, the favored phases from PAW calculations are fully reop¬ 
timized within a force and energy convergence threshold of 10 -5 Ry/bohr and 10 -6 Ry to minimize 
the internal force and energy, respectively. The Troullier-Martins norm-conserving pseudopoten¬ 
tials for Li and C are generated using the FHI98PP code 21 with 2s 1 and 2s 2 2p 2 as valence electrons, 
respectively. The pseudopotentials are then carefully tested by comparing the calculated lattice pa¬ 
rameter and electronic band structure with VASP codes. Convergence test gives the choice of 
kinetic energy cutoffs of 80 Ry, and the Monkhorst-Pack (MP) ^ 22 grids of k-point sampling for 
each favored phases using a grid of spacing 27t x 0.025A L in the Brillouin zone (BZ). A density 
of q mesh with 2n x 0.04A 1 in the first BZ is used in the interpolation of the force constants 
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for the phonon dispersion curve calculations, and a denser mesh of 2n x 0.025A L is used for the 
phonon density of state (PDOS) curve calculations. Subsequently, EPC are calculated in the first 
BZ on the same MP q-point meshes using individual EPC matrices obtained with a grid of spacing 
27t x 0.025A \ All the convergences of the plane-wave basis set and MP sampling are carefully 
examined by employing higher ki netic energy cutoffs and denser grids sets. 

Online Content Supplementary material is available in the online version of the paper. 
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Figure 1 (Color online) Calculated enthalpy difference versus pressure for various 
structures. Enthalpies per CLi 4 unit of structures as functions of pressure with respect to 
the decomposition (C+4Li and (C 2 Li 2 )/2+3Li) enthalpies (black and red lines, respectively), 
which are calculated by adopting the C, Li and C 2 Li 2 structures from references^^and 
at corresponding pressures. 

Figure 2 The structures of CLi 4 under high pressures, a, The P-1 phase at 15 GPa; 
b, P2i phase at 30 GPa; c, C2/m phase at 40 GPa; d, R-3m phase at 60 GPa; e, Fddd 
phase at 150 GPa. Black and blue atoms are C and Li, respectively. 

Figure 3 (Color online) Pressure-induced evolutions of anti-metallization and su¬ 
perconductivity. The evolution from metal to insulator described quantitatively with 
Fermi Surface Filling Parameter (before 70 GPa, pink dash line and olive solid line in 
left part) and Electron Energy Band Gap (after 70 GPa, olive and purple solid lines in 
right part)using ab initio GW calculations, a-e, The 3D FS are shown as upside calcu¬ 
lated by GW method for the corresponding phases at selected pressures, f-g, 3D FS 
of R-3m CLi 4 using ab initio DFT calculation, and its projected FS nesting function £(Q) 
along high-symmetry directions at 60 GPa. The superconducting T c and electron-phonon 
coupling parameter A magnified by ten times as functions of pressure are depicted with 
red and navy solid line in left party, respectively. 

Figure 4 (Color online) Electronic structures of the CLi 4 phases at selected pres¬ 
sures. Blue and black atoms represent Li and C, respectively, a, P2 4 phase at 30G. 
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Left up and down are the PDOS (partial DOS) and integrated PDOS of element Li and C. 
Meddle is 3D-ELF (3 dimensional electron localization function) map with the elf value of 
0.5, and right is the corresponding 2D-ELF (2 dimensional electron localization function) 
slice along the (001) plane, b, Fddd phase at 150G. Left are PDOS and integrated PDOS 
of element Li and C. Meddle represents 3D-ELF with the elf value of 0.5, and right depicts 
corresponding 2D-ELF for the (010) plane. 


21 


